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Abstract. The study of hot plasma expansion in a magnetic field is of interest for many 
astrophysical applications. In order to observe this process in laboratory, an experiment is 
^ ' proposed in which an ultrashort laser pulse produces a high-temperature plasma by irradiation 

J^ of a small target. In this paper an analytical model is proposed for an expanding plasma cloud 

in an external dipole or homogeneous magnetic field. The model is based on the self-similar 
solution of a similar problem which deals with sudden expansion of spheiical plasma into 
^y^ • a vacuum without ambient magnetic field. The expansion characteristics of the plasma and 

^ ' deceleration caused by the magnetic field are examined analytically. The results obtained can 

be used in treating experimental and simulation data, and many phenomena of astrophysical 
and laboratory significance. 
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^^ ] 1. Introduction 

l> 

f^ ' The problem of sudden expansion of hot plasma into a vacuum in the presence of an external 

magnetic field has been intensively studied in the mid-1960s in connection with the high- 
altitude nuclear explosions. It has also been discussed in the analysis of many astrophysical 
and laboratory applications (see, e.g., [1 2| and references therein). Such kind of processes 
KJf ' arise during the dynamics of solar flares and flow of the solar wind around the Earth's 

''T^ , magnetosphere, in active experiments with plasma clouds in space, and in the course of 

C^ ' interpreting a number of astrophysical observations [1-5|. Researches on this problem are 

of considerable interest in connection with the experiments on controlled thermonuclear 
fusion ||6l (a review |[T| summarizes research in this area over the past four decades). 

The expanding plasma pushes magnetic field out which is a consequence of magnetic 
flux conservation. Due to the plasma expansion the surface of the given magnetic flux tube 
inside plasma increases. The conservation of the magnetic flux through the surface of the tube 
leads to a decrease of the local magnetic field and the formation of the diamagnetic cavity. 
Thus even if a magnetic field is nonzero initially inside the plasma it tends to zero during 
plasma expansion. In addition, the plasma is shielded from the penetration of the external 
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field by means of surface currents circulating inside the thin layer on the plasma boundary. 
Ponderomotive forces resulting from interaction of these currents with the magnetic field 
would act on the plasma surface as if there were magnetic pressure applied from outside. 
Thus after some period of accelerated motion, plasma gets decelerated as a result of this 
external force acting inward. The plasma has been considered as a highly conducting media 
with zero magnetic field inside. From the point of view of electrodynamics it is similar to the 
expansion of a superconductor in a magnetic field. An exact self-similar analytic solution for 
a uniformly expanding, highly conducting plasma sphere in an external uniform and constant 
magnetic field has been obtained in |7|. The non-relativistic limit of this theory has been used 
by Raizer 1 8 1 to analyze the energy balance (energy emission and transformation) during the 
plasma expansion. The similar problem has been considered in |J9l within one-dimensional 
geometry for a plasma layer In our recent papers |10| and fTTl we obtained an exact self- 
similar analytic solution for the uniform relativistic expansion of the highly conducting plasma 
sphere or cylinder in the presence of a dipole or homogeneous magnetic field, respectively. 

As mentioned above the previous treatments 17-11] were obtained assuming somewhat 
idealized situation: uniform expansion, zero magnetic field and thermal pressure inside 
plasma, etc. Obviously such simplified models are not capable to describe correctly the 
plasma expansion dynamics where an essential deceleration (or acceleration) of the plasma 
boundary may occur More realistic models for the plasma self-similar expansion have 
been developed, for instance, in lfT2l - [T4l (see also references therein) for planar (one- 
dimensional) 1 12 1 and cylindrical fl3"141 expansions employing ideal magnetohydrodynamic 
(MHD) equations. It should be noted that the ideal MHD is well justified since the typical 
experimental parameters are such that the expanding plasma is coUisionless and the dissipative 
effects are negligible HJ-Hl. 

In the present paper we study the expansion of the spherical plasma cloud in the presence 
of a dipole or homogeneous magnetic field taking into account the thermal effects. We employ 
the self-similar solution of a similar problem which deals with expansion of spherical plasma 
into a vacuum without magnetic field. Unlike the previous studies llT2Hl4l we concentrate 
here on the dynamics (deceleration and/or acceleration) of the sharp plasma boundary. We 
found an analytical solution which can be used in analyzing the recent experimental and 
simulation data (see, e.g., |[T]|2l and references therein). 

2. Theoretical model 

Usually the motion of the expanding plasma boundary is approximated as the motion with 
constant velocity (uniform expansion). In the present study a quantitative analysis of plasma 
dynamics is developed on the basis of one-dimensional spherical radial model. Within the 
scope of this analysis the initial stage of plasma acceleration, later stage of deceleration and 
the process of stopping at the point of maximum expansion are examined. 

We consider a coUisionless magnetized plasma expanding into vacuum. The relevant 
equations governing the expansion are those of ideal MHD ITSl assuming that the 
characteristic length scales for the plasma flow are much larger than the Debye length and 
the Larmor radius of the ions. Thus 



^ + (v.V)v 



P 



— [V X H] X H - VP, 



^^ = V X [v X H] , V • H = 0, (1) 
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where p, v, H, and P are the mass density, the velocity, the magnetic field, and the pressure 
respectively. These equations must be accompanied by the adiabatic equation of state P ^ p^ 
(7 = Cp/Cy > 1 is the adiabatic index, Cp and Cy are the heat capacities at constant 
pressure and constant volume, respectively) and the equation of conservation of entropy, 
which expresses the fact that the plasma dynamics is adiabatic in the absence of dissipation. 
Using the thermodynamic relation between entropy, pressure, and internal energy as well as 
equation ([TJ the equation for pressure reads 1 16] 

dP 

— + (v-V)P + 7P(V-v)-0. (2) 

The equation of conservation of energy is derived from the system of equations ^ and 
is given by ifTSl 

de 

^ + V-Q-0, (3) 

where e and Q are the energy and energy flux densities, respectively, with 
_piP_ P IP 

Q = v(^ + -2^)+-l[Hx[vxH]]. (4) 

Here the energy e consists of the kinetic (first term), the internal (second term) and the 
magnetic field (third term) energies. The last term in the energy flux density Q represents the 
Poynting vector (let us recall that in ideal MHD the electric field is given by E = — i [v x H]). 
With the theoretical basis presented so far, we now take up the main topic of this paper 
This is to study the dynamics of the plasma boundary expanding in an ambient magnetic field. 
Consider the magnetic dipole p and a plasma spherical cloud with radius a{t) located at the 
origin of the coordinate system. The dipole is placed in the position rg from the center of 
the plasma cloud (a(i) < tq). The orientation of the dipole is given by the angle 9p between 
the vectors p and Tq. We denote the strength of the magnetic field of the dipole by Ho(r). 
The energy, which is transferred from plasma to electromagnetic field is the mechanical work 
performed by the plasma on the external magnetic pressure i^Q (r)/87r. Taking into account 
this effect and the energy conservation ^ integrated over the spherical plasma volume flp 
(with ^ r ^ a(t)) the equation of balance of plasma energy is as follows: 

/ P{r,t)r^dr + 2tt / p{r,t)v^{r,t)r'^dr (5) 

7 - 1 Jo Jo 

where il is the volume of the spherical shell aq ^ r ^ a{t), Oq = a(0) (a{t) ^ oq) and 
Wq are the initial radius and energy of the plasma. When the plasma cloud is introduced 
into a background magnetic field, the plasma expands and excludes the background magnetic 
field to form a magnetic cavity. The magnetic energy of the dipole in the excluded volume 
is represented by the last term in equation (|5]). Initial plasma velocity is supposed to be 
v{r, 0) — v„i{r/ao) at r ^ ao and v{r, 0) = at r > qq, where u,„ is the initial velocity of 
the plasma boundary (v„i = a(0)). 

The obtained energy balance equation can be effectively used if profiles of velocity 
v{r, t), pressure P{r, t), and mass density p{r, t) are known functions of the plasma radius 
a{t). We will take these dependencies from the solution of a similar problem which deals with 
sudden expansion of spherical plasma into a vacuum without ambient magnetic field I1I6J . 
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The simplest class of solutions available in this case are so-called self-similar solutions. 
They are realized under the specified initial conditions. We will set the initial conditions 
with a parabolic distribution of pressure and mass density, which describe hot and dense 
initial plasma state with sharp boundary localized at r — oq. The self-similar solutions are 
characterized by a velocity distribution linearly dependent on r (see, e.g., llT6l ). At r ^ a{t) 

v{r,t)^r^ (6) 

a{t) 

where unknown a{t) is the radius of sharp plasma boundary while a{t) is the velocity of the 
boundary. The specification of the mass density profile at r ^ a{t) is given by 



P {r, t) - 



1 



m 



(7) 



7r3/2r (1 + q) a3 [t) 

and equation (|6]l for the velocity, automatically satisfies the continuity equation in ([T]i for an 
arbitrary function a{t) and for an arbitrary parameter q. Here M ~ const is the total mass of 
plasma cloud and r(z) is the Euler function. Substitution of w(r, t) into the entropy equation 
(|2| gives at r ^ a{t) the following solution for the pressure 



P(r,i)=p„ 



ao 



a{t) 



13')' r y,2 
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{t) 



(8) 



where s is an arbitrary parameter and Pmax is the thermal pressure at the center of the spherical 
plasma cloud at t = 0. In addition the quantities v{r^ t), p{r, t) and P{r, t) vanish at r > a{t), 
v{r, t) = p(r, t) ~ P{r, t) — 0. Substituting ^ and (O into the fluid equation of motion ([T]) 
and neglecting the term involving the magnetic field H yields a second-order differential 
equation governing the motion of the plasma boundary a{t). The problem considered is 
not isentropic in general except the case when s — q-f. In the latter case of the isentropic 
expansion the equation of state is given by Pp '' — const. Throughout in this paper we will 
assume that q ^ and s ^ 0. The self-similar solution given by equations ©-([HI) has been 
considered by many authors. We refer to |fT6l (see also the references therein) for details and 
for some applications of this solution. 

There are several limitations and shortcomings of the self-similar model considered 
above. It should be emphasized that in general the plasma expansion process with and without 
ambient magnetic field is not self-similar since there is a characteristic length scale, that is, 
the initial radius of the plasma qq. However, during the later stage (when a{t) ^ ao) the 
plasma expansion asymptotically approaches to the self-similar regime since the role of the 
radius aq becomes less important. This property is supported by the asymptotic solution of 
the system of partial differential equations ([T]i and Q lfT6HT7l as well as by the numerical 
simulations [16] for the free expansion (i.e. at H = 0). In addition, it has been shown that 
at the beginning stage of the expansion the deviation of the solution ©-([H]) from non-similar 
one is small and is additionally reduced due to the spatial average in the equation of energy 
balance Q. 

Equations ©-([Hi are an exact solution in the case of expansion into a vacuum without 
magnetic field. However, equation (O does not satisfy the boundary condition, P{a{t), t) = 
Hq/Stt, which is imposed in the case of expansion into an ambient magnetic field. Therefore, 
in the present case with nonzero magnetic field for the validity of the equations (|6]l-([8j most 
critical is the domain close to the boundary of the plasma, r ~ a{t). On the other hand 
if the magnetic pressure is smaller than the plasma average pressure P, Pmag ^ P, the 
difference between the exact solution in the magnetic field and free expansion model is 
small and is localized in a narrow area near the surface of the cloud. These deviations are 



Self-similar analytical model of the plasma expansion in a magnetic field 5 

additionally reduced due to the integration in the equation of energy balance. Estimating 
accuracy of the free expansion model, one should take into account that the long stage of 
plasma deceleration corresponds to a high expansion ratio, a{t)/aQ ^ 1. Average plasma 
pressure drops significantly and plays no role in energy balance equation (|5]l during this stage. 
In accordance with the above boundary condition local pressure near the plasma edge must 
be equal to the magnetic pressure outside. It causes deviation from the profile equation (|7]i 
and accumulation of plasma in this area. This is confirmed independently by the numerical 
simulations ITSl . In the limiting case when all plasma is localized near the front, one can 
expect an increase of the kinetic energy and longer stage of plasma deceleration as compared 
with the self-similar expansion model. 

Another critical domain for the violation of the self-similar solution ©-(ISll is the final 
stage of expansion when the plasma is fully stopped by the magnetic field. As mentioned 
above the average plasma pressure is strongly reduced compared to the magnetic pressure and 
the equations (|6]l-(l8]l clearly become invalid in this case. However, if the critical time interval 
At, where P < Pmag, is much smaller than the typical time scale of the plasma flow (up to 
the full stop) the contribution of this interval to the overall plasma dynamics is negligible and 
the use of the self-similar solution is justified. 

In the case of dipole magnetic field the volume integral in the last term of equation (|5]i 
has been evaluated in [J0|. The result reads 

^^dr = ^[Qir^x{t))-Q{7^)], (9) 

where rj = ap/ro < 1, x{t) = a{t)/ao (note that a{t)/ro = rix{t) < 1), and 

Q{V)= ,^ \. [v{l-v'){3cos'0p-l) (10) 

„ Q ., 2 /, M 3cos^6'„ - 1 , 1 + ?7 

+8i]^ (1 + cos^ Op)] -^ In '-. 

^ '■' 2 1 ^ ^ 

Substituting equations ©-([S]) into ^ and integrating over r yields first-order differential 
equation for a{t), which already satisfies initial condition a(0) = Vm, 

i\r) + -^+a[Q{f^x{T))^Q{rj)]^l. (11) 

Here two dimensionless quantities are introduced 



which determine the magnetic and the thermal energies, respectively, in terms of the total 
initial energy Wg. The latter is easily obtained from equation (|5]) and reads 

3Mvl , TT^/ V,,,a3 r(l + .) 

^° = 2(5T^+ 7-1 r(| + .)- ^'^^ 

From equations ( fT2] i and (fT3] i it is seen that /3 < 1. New dimensionless variables are 
introduced as follows: x{t) = a{t)/aQ, r = t/to, to — ao/um, where Um ~ [2(5 + 
2g)Wo/3M]^/^ is the velocity of plasma expansion, achieved asymptotically at t -^ cx) in the 
case of expansion into a vacuum without magnetic field (i.e. at a = 0). 

The total energy of the plasma cloud at time t is obtained from equation (fTTb 

2 

Wit) = Wo^^[Q{nx{t))^Q{r^)]. (14) 
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Note that the function Q(?/) monotonically increases with the argument and the plasma cloud 
energy decreases with time. 

Consider also the case of uniform magnetic field when Hq — const. In this case the 
volume integral in equation (|9]l is replaced by {Hq/Q) [a^ (i) — Og ) and the differential equation 
(fTTT l for the plasma boundary reads 

*'M + ^3(^ + ^[■^'(^) -!] = !' (15) 

where a = W^ag/Wo, W,„ag = (47rao/3)Pi„ag is the initial magnetic energy in the plasma 
volume, and Pmag = Hq/8tt is the magnetic field pressure. 

Equations (fTTT l and ( fTSI l coincide with the equation of the one-dimensional motion of 
the point-like particle in the potential U{x) which is determined by second and third terms of 
equations ( fTTT l and (fTsT l. The distance Xg of the plasma cloud motion up to the full stop (the 
stopping length) at the turning point is determined by U{xs) = 1. In particular, it is easier to 
obtain the stopping length in the case of homogeneous magnetic field and at vanishing thermal 
pressure {(3 — 0). Then from equation (fTST i one obtains the equation of motion 

x^ = l + a-ax^. (16) 

It is seen that in this case the stopping length is given by Xs = (1 + l/aY^^. The solution of 
equation (fTST i can be represented in the form 

where J^{z) = F (i i; |; z) and the latter is the hypergeometric function. Substituting in 
equation ( fTTj ) x{t) — Xg we obtain the corresponding stopping time as a function of the 
magnetic field and the plasma kinetic energy 



i = 



4«^(^V^fi 



(17) 



t '' 



Vo-+l 



1/3 



1 



(18) 



Here C = J'(l) = y^T (|) /T (|) ~ 1.4 is a constant. At vanishing (a < 1) 
and very strong (a ^ 1) magnetic fields the stopping time becomes tg ~ CRm/vm — 
C{ao/vm)o'^^^'^, ts ~ 2Rf^/3vmaQ ~ (2ao/3wm)c^^, respectively, where the radius 
Rm — (QWo/H^y^'^ is obtained by equating the initial kinetic energy Wo of an initially 
spherical plasma cloud to the energy of the magnetic field that it pushes out in expanding to 
the radius Rm- It is worth mentioning that in the case of weak magnetic field, cr <C 1, and at 
vanishing thermal pressure (/? = 0) the stopping time does not depend on the initial plasma 
radius, t, - (Af/w™Pmag)'/^. 

We now turn to the general equations determined by (fTTT l and ( fTsT l. At the initial stage of 
plasma expansion (i ^ Iq) from these equations we obtain 

x{t)c.l+^-^ + h(l-)\ (19) 

flo 4 \toJ 

where h — f3{'-f — 1) — k, k — %riQ'{ri) and n — a for the dipole and homogeneous magnetic 
fields, respectively. Here the prime indicates the derivative with respect to the argument. Thus 
at the initial stage the plasma cloud may get accelerated or decelerated depending on the sign 
of the quantity h (in other words on the relation between thermal and magnetic pressures). 
For instance, in the homogeneous magnetic field the acceleration occurs when Pmax > Pc, 
where 



4 r( 



5 



S 



Pc — ^ I— -n fT , \ -'mag (^U) 

6y/TT r (1 + S) 
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Figure 1. The dynamics of the plasma cloud expanding in a dipole magnetic field. Shown ai'e 
the scaled radius a{t)/ao (left panel) and the velocity d{t)/um (right panel) of the plasma 
boundary vs time (in units of to) at 7 = 5/3, = 0.5, t; = 10^^, 9p = and for « = 0.1 
(solid lines), a = 100 (dashed lines) and a = 2 X 10* (dotted lines). 



(i.e. at h > 0) and continues until x{t) reaches some value Xc > I given by Xc = 
(Pmax/Pc)^^'^^- The time interval ^ t < tc of the acceleration is determined from the 
equation of motion ( fTSl ). The critical radius Xc and time tc correspond to the beginning of 
plasma deceleration. Further plasma motion at t > tc is an expansion with slowing-down 
velocity. It ends up at the turning point which corresponds to the maximum of expansion, 
U{xs) = 1. However, in the opposite case of the low thermal pressure with pmax < Pc the 
plasma systematically get decelerated in the whole time interval of its dynamics. 

A characteristic stopping time of plasma motion up to the full stop at the turning point is 
given by the integral of the equations ( fTTl i and ( fTSl l 



dy 



1 



U'ixs) 



(21) 



\/l - U{y) 

Calculating time tg needed for plasma to reach this point, one can simplify the integrand 
taking into account that the main contribution comes from the vicinity of upper limit of 
integration. This approximation is expressed by the second part of equation (ISTT l. In 
the case of weak and homogeneous magnetic field this yields universal expressions, tg ^ 



{M/umPm^^f^ and , 



ji^. It is worth mentioning that in the case of weak magnetic 



field the stopping time and length do not depend on the initial plasma radius but depend on 
the thermal pressure (or temperature) (cf. these relations with those obtained above). At very 
strong magnetic fields, a^ ~ oq -I- (l/2)umis and tg ^ Mvm/a^Pma.g, and the stopping 
characteristics of the plasma essentially depend on the initial radius but are now independent 
on the thermal pressure. The similar estimates can be found for the dipole magnetic field. 
However, we note that the latter case significantly differs from the homogeneous field situation 
considered above. Since in the vicinity of the dipole the magnetic field is arbitrary large the 
stopping length cannot naturally exceed ro for any thermal energy of the plasma {x{t) < l/rj 
in (fTTTl). For a weak magnetic field this simply yields Og ~ ro and tg ~ ro/um. 

As an example in figure[T]we show the results of model calculations for the normalized 
radius a{t)/aQ (left panel) and the velocity a{t)/u,-a (right panel) as a function of time (in 
units of io) at 7 = 5/3, /3 — 0.5, rj — 10~^, 6p = Q and for different values of the parameter 
a. In this figure the dimensionless strengths a of the dipole magnetic field are chosen such 
that the coefficient h in equation ( fT9] l is positive, /i > 0, for solid and dashed lines and 
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Figure 2. The normalized stopping length as(ci)/aa (left panel) and the stopping time 
ts{ct)/to (right panel) of the plasma cloud expanding in a dipole magnetic field vs the 
normalized dipole magnetic field a at 7 = 5/3, 7] = 10~^, 6p = and for /3 = 0.1 
(solid lines), f} = 0.5 (dashed lines) and (5 = 0.9 (dotted lines). 



negative, h < 0, for dotted lines. From the right panel of figure [T] it is seen that at h > 
(solid and dashed lines) there is a short initial period of acceleration, ^ t < to, when the 
plasma boundary is accelerated according to equation ([% . During this period (which is only 
weakly sensitive to the magnetic field strength) the dimensionless radius a{t)/ao increases up 
to 2 — 3, and at io ^ ^ < ic almost all initial total energy Wq is transferred into kinetic energy 
of free radial expansion at constant velocity ^ u„i. As expected (see above) the time tc is 
reduced with increasing the strength of the magnetic field and the free expansion period is 
shorter for larger a. The further increasing the strength of the magnetic field (figure[Tl dotted 
line) results in a plasma dynamics with systematically slowing-down velocity. 

For the same set of the parameters 7, 77 and dp in figure |2] it is shown the normalized 
stopping length (left panel) and the stopping time (right panel) of the plasma cloud as a 
function of the dimensionless strength a of the dipole magnetic field for some values of the 
normalized plasma thermal pressure (3. It is seen that the stopping length and time decrease 
with the strength of the magnetic field and practically are not sensitive to the variation of the 
plasma thermal pressure. 

Note that at otherwise unchanged parameters the strength of the dipole magnetic field is 
maximal at the orientation 9p = and monotonically decreases with dp. For instance, the 
strength Hq{0) of the dipole magnetic field at the center of the plasma cloud is reduced by a 
factor of 2 by varying the dipole orientation from 6p — to dp = 7r/2. Therefore the effect 
of the magnetic field shown in figures [T] and |2]is weakened at the orientation 9p — tt/2 of the 
dipole. In particular, this results in a larger stopping lengths and times than those shown in 
figure |2] 

In this paper we have assumed that the self-similar solutions ©-dSJ and hence the 
shape of the plasma remain isotropic during plasma expansion. Such assumption is not 
evident since the distribution of the magnetic pressure on the plasma surface is anisotropic 
in general. Thus one can expect strong deformation of the initially spherical plasma surface. 
Consider briefly this effect assuming constant magnetic field Hq. Taking into account the 
effect of the induced magnetic field the total magnetic pressure on the plasma surface is 
Pm{&) = (9/4)Pniag sin^ 9, where 9 is the angle between the radius vector r and Hq fS'TOl. 
The magnetic pressure vanishes at 6* — 0, vr and is maximal at 6* = tt/2. Therefore, 
one can expect that the spherical shape of the plasma is deformed into the ellipsoidal one 
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which is elongated along the magnetic field lines. This effect is supported by the numerical 
simulations iHHSl. Thus, the isotropic self-similar model considered here is valid as long as 
the deformation of the plasma shape is much smaller than the typical length scale for the 
plasma flow. However, while the solution ©-([HI) is violated in the case of strong anisotropy 
the plasma dynamics can be described (at least qualitatively) by equations (fT5])-(l2T]) replacing 
the isotropic magnetic pressure Pmag by the angular dependent one Pm{9). The expansion 
radius a{t, 0) as well as the stopping length asiO) and time t8{d) should be also anisotropic 
in this case. 

3. Conclusions 

An analytical self-similar solution of the radial expansion of a spherical plasma cloud in the 
presence of a dipole or homogeneous magnetic field has been obtained. The analysis of the 
plasma expansion into ambient magnetic field shows that there are processes of acceleration, 
retardation and stopping at the point of maximum expansion that are very distinct and 
separated in space and time. The scaling laws obtained are, in general, the functions of two 
dimensionless parameters, a (or a for constant magnetic field) and /?, which can be varied by 
means of the choice of the external magnetic field, the thermal pressure and the initial energy 
of the plasma. It allows to test the different regimes of the plasma dynamics in a wide range 
of external conditions. 

We expect our theoretical findings to be useful in experimental investigations as well 
as in numerical simulations of the plasma expansion into an ambient magnetic field (either 
uniform or nonuniform). One of the improvements of our present model will be the derivation 
of the dynamical equation for the plasma surface deformation. In this case it is evident that 
the problem is not isotropic with respect to the center of the plasma cloud (r = 0) and a 
full three-dimensional analysis is required. A study of this and other aspects will be reported 
elsewhere. 
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